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Abstract 

In this paper we present an approach towards the comprehensive anal- 
ysis of the non-integrability of differential equations in the form x = f(x, t) 
which is analogous to Hamiltonian systems with 1 + 1/2 degree of free- 
dom. In particular, we analyze the non-integrability of some important 
families of differential equations such as Painleve II, Sitnikov and Hill- 
Schrodinger equation. We emphasize in Painleve II, showing its non- 
integrability through three different Hamiltonian systems, and also in Sit- 
nikov in which two different version including numerical results are shown. 
The main tool to study the non-integrability of these kind of Hamiltonian 
systems is Morales-Ramis theory. This paper is a very slight improvement 
of the talk with the same title delivered by the author in SIAM Conference 
on Applications of Dynamical Systems 2007. 

Key words and phrases. Hill-Schrodinger equation, Morales-Ramis 
theory, non-autonomous Hamiltonian systems, non-integrability of Hamil- 
tonian systems, Painleve II equation, Sitnikov problem, virtually abelian 
groups. 
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1 Introduction 

In this section we present the necessary theoretical background to understand 
the rest of the paper. 
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1.1 Differential Galois Theory 

Our theoretical framework consists of a well-established crossroads of Dynamical 
Systems theory, Algebraic Geometry and Differential Algebra. See [20] or [34] 
for further information and details. Given a linear differential system with 
coefficients in C(t), 

z = A(t)z, (1) 

a differential field L D C(t) exists, unique up to (C(t)-isomorphism, which con- 
tains all entries of a fundamental matrix if? = . . . ,i/>„] of {2}. Moreover, 
the group of automorphisms of this field extension, called the differential Galois 
group of ([T]), is an algebraic group G acting over the C- vector space (if) l , . . . , tj) n ) 
of solutions of {T]) and containing the monodromy group of ([I]) . 

It is worth recalling that the integrability of a linear system ([I]) is equivalent 
to the solvability of the identity component G° of the differential Galois group 
G of {T]) - in other words, equivalent to the virtual solvability of G. 

It is well established (e.g. [H]) that any linear differential equation system 
with coefficients in a differential field K 

± ( \ _ ( o(t) b(t) \(& 

dt\i2 \ C(t) d(t) U 



(2) 

by means of an elimination process, is equivalent to the second-order equation 




bit)c(t) a(t)d(t) j C = o, ( 3 ) 

where £:=£].. Furthermore, any equation of the form z — 2pz — qz = 0, can 
be transformed, through the change of variables z = yef p , into y = —ry, r 
satisfying the Riccati equation p = r + q + p 2 . This change is useful since it 
restricts the study of the Galois group of y = —ry to that of the algebraic 
subgroups of SL2 (C). 

A natural question which now arises is to determine what happens if the 
coefficients of the differential equation are not all rational. A new method was 
developed in [3] , in order to transform a linear differential equation of the form 
x = r(t)x, with transcendental or algebraic non-rational coefficients, into its 
algebraic form - that is, into a differential equation with rational coefficients. 
This is called the algebrization method and is based on the concept of Hamil- 
tonian change of variables [3] . Such a change is derived from the solution of a 
one-degree-of-freedom classical Hamiltonian. 

Definition 1.1 (Hamiltonian change of variables) A change of variables 
t = t (t) is called Hamiltonian if (r(t),f(t)) is a solution curve of the au- 
tonomous Hamiltonian system Xh with Hamiltonian function 

p 2 ~ 

H = H{r,p) = y + V(t), for some V € C(r). 
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Theorem 1.2 (Acosta-Blazquez algebrization method |3j) Equation x = 
r(t)x is algebrizable by means of a Hamiltonian change of variables r = r(t) if, 
and only if there exist /, a such that -4- (In a) , ^ S C(r), where 

f{r{t))=r{t), a{r)^2{H~V{r)) = (Tf. 
Furthermore, the algebraic form of x = r(t)x is 

Id, \di /f\ 

o3Z ln « IT" - * = 0. □ (4) 



The next intended step, once a differential equation has been algebrized, is 
studying its Galois group and, as a causal consequence, its integrability. Con- 
cerning the latter, and in virtue of the invariance of the identity component 
of the Galois group by finite branched coverings of the independent variable 
(Morales- Ruiz and Ramis, [331 Theorem 5]), it was proven in [3J Proposition 1] 
that the identity component of the Galois group is preserved in the algebrization 
mechanism. 

The final step is analyzing the behavior of t = oo (or t = oo) by studying 
the behavior of n = through the change of variables r\ = 1/t (or n = 1/r) 
in the transformed differential equation, i.e. t = oo (or r = oo) is an ordinary 
point (resp. a regular singular point, an irregular singular point) of the original 
differential equation if, and only if , rj = is one such point for the transformed 
differential equation. 



1.2 Morales-Ramis Theory 

Everything is considered in the complex analytical setting from now on. The 
heuristics of the titular theory rest on the following general principle: if we 
assume system 

y = X(z) (5) 

"integrable" in some reasonable sense, then the corresponding variational equa- 
tions along any integral curve r = {z (t) : t E 1} of © , defined in the usual 
manner 

€ = (VE r ) 

must be also integrable - in the Galoisian sense of the last paragraph in ll.ll We 
assume T, a Riemann surface, may be locally parametrized in a disc / of the 
complex plane; we may now complete T to a new Riemann surface T, as detailed 
in [UJ §2.1] (see also [501 §2-3]), by adding equilibrium points, singularities of 
the vector field and possible points at infinity. 

The aforementioned "reasonable" sense in which to define integrability if 
system ([5]) is Hamiltonian is obviously the one given by the Liouville-Arnold 
Theorem, and thus the above general principle does have an implementation: 

Theorem 1.3 (J. Morales-Ruiz & J. -P. Ramis, 2001) Let H be an n-degree- 
of-freedom Hamiltonian having n independent first integrals in pairwise involu- 
tion, defined on a neighborhood of an integral curve T. Then, the identity com- 
ponent Gal (VEp)° is an abelian group (i.e. Gal (VEp) is virtually abelianj. 
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See [331 Corollary 8] or [30] Theorem 4.1] for a precise statement and a proof. 

1.3 Non Autonomous Hamiltonian Systems 

Non-autonomous Hamiltonian systems (NAHS) on symplectic manifolds have 
long been the subject of study, and appear in a most natural way in Classical 
Mechanics and Control Theory, e.g. pQ, [5], [17], [19], [18], [26], pH], [32]. 
We consider NAHS of the form 

2 

H = H(q 1 ,p 1 ,t) = ^-+V(q 1 ,p 1 ,t), (6) 

is an NAHS with 1 + 1/2 degree of freedom. It is well-known (e.g. that 
(JHJ) can be included as a subsystem of the Hamiltonian system with two degree 
of freedom given by 

p 2 

H = H(qi,q 2 ,pi,p 2 ) = — + V(qi,pi,q 2 ) + p 2 , (7) 

where q 2 andp2 are conjugated variables, i.e. p 2 = —H + k, where k is constant, 
and q 2 — t. Furthermore, p 2 is easily seen to be the offset or counterbalancing 
energy of the system ([2B], [2"8]). 

Also worth mentioning are some recent results on canonical transformations 
in the extended phase space [28] , [29] , [30] , [33] , as well as on definitions and con- 
sequences of "integrability" under such circumstances or generalizations thereof, 
even for non-Hamiltonian systems ([10]. [9], [12], [16]) which we will not delve 
into further at this point. 

2 Main results 

Consider the differential equation 

x = f(x,t), (8) 

with particular solution x = x(t). We will henceforth order our choice of posi- 
tions as qi = x and q 2 = t, thus yielding a Hamiltonian system given by 

H = — -F{q 1 ,q 2 ) , F qi (qi,q 2 ) = ^ = f(qi,q2)- 

2 oqi 

Equation ((5J) is obviously equivalent to Hamilton's equations for H, 
q 1 =p 1 =H pi pi = -H qi = f(qi,q 2 ); 

this NAHS is included as a subsystem of linked to H := H + p 2 , such as 
in equation ([7]). Assuming x(t) = q\(i) to be a solution of {5J and 52 (i) = we 
obtain an integral curve r = {z (<)} of H, where 

z{t) := {qi{t),q2{t), Pl (t),p 2 (t)) = {q x {t),t,qi{t),-H{t)) . 

We may now introduce our first main result: 
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Theorem 2.1 Let T be an integral curve of X^ such as the one introduced 
above. If X^ is integrable by means of rational or meromorphic first integrals, 
then the Galois group of 

£ =(/«*(«!, fcOIrK, (9) 

is virtually abelian. 

Proof. The Hamiltonian field Xq is given by = (pi, 1, f{q\, (72), F q2 (qi, q2]) T 
The variational equation VEr along T is 





/a n 




/ 1 \ 
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~dt 






/gi (91,92) /g 2 (91,92) 
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More precisely, equation (fT0|) is 
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(10) 



(11) 



Hence £2 = k, where k is constant. Assuming k = 0, the normal variational 
equations (NVE r , see §1], [2J §4.3], [201 §4.1.3]) for H are given by 



dt 



6 



1 

/ 9l (9i,92)|r 



6 



(12) 



and solving (fT2)) we can obtain £4- System (TTJJ is equivalent to equation (JOJ, 
where £ = £1. In virtue of [201 Proposition 4.2], the virtual commutativity of 
Gal (VEr) implies that of Gal (NVEr); this, coupled with Theorcm ll.31 implies 
that if Xfi is rationally or meromorphically integrable, then the Galois group 
of the equation (J0J is virtually abelian. ■ 



Remark 1 The above disjunctive between meromorphic and rational Hamilto- 
nian integrability is related to the status of t = 00 as a singularity for the normal 
variational equations. More specifically, and besides the non-abelian character 
of the identity component of the Galois group, in order to obtain Galoisian ob- 
structions to the meromorphic integrability of H the point at infinity must be a 
regular singular point of ([9j. On the other hand, for there to be an obstruction 
to complete sets of rational first integrals, t — 00 must be a irregular singular 
point. 

Corollary 2.2 If the Galois group of the differential equation £ — fc(i)£, is not 

2 2 ^ 

virtually abelian, then, defining H := Q — fc(9a) 2 an< ^ H := H + P2, X^ is not 
integrable by means of meromorphic or rational first integrals. □ 
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Remark 2 If the Galois group of the equation £ = k(t)£ is the Borel group 
G = C* k C (hence connected, solvable and non-abelian), is neither mero- 
morphically nor rationally integrable, although it is still possible to solve the 
equation - as well as, ostensibly, the NAHS Xh- 

In expectance of the following Corollary consider, for any g(x), a(t) and 
a — ao given, the equation 

x = g x (x)(g(x) + a(t)) + a, a € C, (13) 
having a certain known particular solution q\ = q\ (t). Let 

„ pi (g(gi) + «(g2)) 2 

H= — aqi, q 2 = t 

be a Hamiltonian linked to (|13j) . and H := H + p2 its autonomous completion. 

Corollary 2.3 If X^ is integrable through rational or meromorphic first in- 
tegrals then, along the integral curve T = {z(t) — (qi(t),t,qi(t),—H(t))}, the 
Galois group of the equation 

i= (9 2 gi {Q^+9niM(9{qi)+a(q2)))\ r ^ (14) 

is virtually abelian. □ 

Now, keeping the above hypotheses for g{x), a(t), a = cxq, equation (fl3|) and 
x = qi(t), let us define the Hamiltonian system 

„2 

H = H + p 2 , H = ^-(g(q 1 )+a{q 2 ))p 1 -(a + a q2 {q 2 ))q 1 , a £ C. (15) 

It only takes the following simple calculation to prove that Hamiltonian H is 
indeed linked to (|13[) in the manner expected, i.e. as introduced in the beginning 
of Section [2] We have 

x = g x (x)(g(x) + a(t)) + a 

= g x {x)x + g x (x)g(x) + g x {x)a(t) + a + a(t) - g x (x)x - a(t) 
= g x (x)y + a + a(t)-j i (g(x)+a(t)), 

where y = x + g(x) + a(t), and requiring x and y to be conjugate variables 
implies the following system, equivalent to (fT3|) . 

£ = y - (g(x) + a(t)) = h v 

y = g x {x)y + a + a(t) = -H x . 
A straightforward parallel integration, 

y 2 

H = — - (g[x) + a(t)) y + hi(x,t) = h 2 (y,t) - g(x)y - (a + a(t))x, 
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along with the definition of 



hi(x,t) = -(ot + a(t))x, h 2 {y,t) = 



a(t)y, 



as well as (51, <fe>Pi) = (x,t,y), yields the autonomous Hamiltonian system 
introduced in IT5l). 



Theorem 2.4 If Xn is integrable by means of rational or meromorphic first 
integrals, then, along T = {z (t) = (qi(t),t,pi(t), — H(t))}, the Galois group of 
the equation 

£= (dlAll) +P-L9qiqi{qi) - 9q 2 qi(<ll))\ T t, (16) 

is virtually abelian. 

Proof. We may proceed as in the proof of Theorem 12.11 The Hamiltonian 
field Xg is given by 



Pi 



(g(qi) +a(q 2 )) \ 
1 

9qi{qi)Pl + (#2)) 

V a qi(qi)P\ +a q2q2 {q 2 )qi ) 



The variational equation VEr along T = {(qi(t),t,pi(t), —H(t))} is 



1 = 



/ 


99(91) 


da(q 2 ) 


dqi 


dqi 












d 2 g(qi) 


d 2 a(q 2 ) 




dqi 
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d 2 a(q 2 ) 
dqi 


I 2 a(q 2 ) 



d 3 a(q2) 



1 


\ 








99(91) 





dqi 


da(q 2 ) 
dqi 


J 



(17) 



where £ = (£1, £2, £3, ^4) T ■ £,2 = k G C; in particular, k = renders system (fT7|) 
equal to the following: 



-5gi(9l)l r 6 + 
0, 

PiSW(?i)l r £i + -991(91)^^3, 

a 92 92 



(18) 



NVEr corresponding to 



d 
df 



6 

6 



-991(91) 1 

Pl9qiqi(qi) 9qMl) 



6 
6 



(19) 



and solving IT^l) we obtain £4. Using equations and ([3]), system (fig]) is equiv- 
alent to equation (TT6]) , where £ = £1. Again in virtue of [201 Proposition 4.2] as 
in Theorem 12. li the integrability of Hamiltonian H by means of meromorphic 
or rational first integrals implies the virtual commutativity, of the Galois group 
of equation (fTH)) . ■ 
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Remark 3 We observe that if x = f(x, t) with particular solution x = x(t) is 
the same for Theorem 12.11 Corollary 12.31 and Theorem 12. 4i then we obtain the 
same NVEp (equations (fT4")l and (jTHJ) are equivalent), despite the fact that 
their respective linked Hamiltonian systems can be expressed differently. 

3 Examples 

In this Section, and in application of Theorems 12.11 and 12.41 as well as of Corol- 
laries 12.21 and 12. 3[ we analyze the non-integrability of the Hamiltonian systems 
corresponding to the following differential equations: 

1. Hill-Schrodinger equation: x = k(t)x, 

2. Painleve II equation: 'x — 2x 3 + tx + a, 

3. The differential equation: x — — \ + ct, 

4. The Sitnikov problem: x = _ii^££££% r ( t ) = i-ecoBt _ 

In order to analyze normal variational equations, a standard procedure is us- 
ing Maple, and especially commands dsolve and kovacicsols. Whenever the 
command kovacicsols yields an output "[ ]", it means that the second-order 
linear differential equation being considered has no Liouvillian solutions, and 
thus its Galois group is virtually non-solvable. For equations of the form y = ry 
with r € <C(x) the only virtually non-solvable group is SL 2 (C). In some cases, 
moreover, dsolve makes it possible to obtain the solutions in terms of special 
functions such as Airy functions, Bessel functions and hypergeometric func- 
tions, among others (0). There is a number of second-order linear equations 
whose coefficients are not rational, and whose solutions Maple cannot find by 
means of the commands dsolve and kovacicsols alone; this problem, in some 
cases, can be solved by a previous algebrization procedure. 

3.1 Hill-Schrodinger equation x = k(t)x 

This example corresponds to Corollary 12. 21 For k > 0, e ^> 0, V n & polynomial 
of degree n with "P„(0) = 1 and k(t) given by 

k(t) = /ce~ £t , k(t) = kV n (et), 

k(t) = fe(l + sinh(ei)), k(t) = k(l + sin(et)), 

k(t) = k(l + cosh(et)), k(t) = fc(l + cos(ef)). 

Xfi is non-integrable by means of rational first integrals. 

The integrability of equation 'x — k(t)x for these examples has been deeply 
analyzed in [3], 
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3.2 The Sitnikov problem 

The Sitnikov problem is a symmetrically configured restricted three-body prob- 
lem in which two primaries with equal masses move in ellipses of eccentricity e 
in a plane tt\ , and an infinitesimal point mass moves along the line tt ± . See [6] , 
[15] , [2"3"] . [2"5] . [55] . [55] . [T3] for more details. The motion of the infinitesimal 
point mass is given by the following differential equation 

z + m = 0, (20) 

(r 2 (t) + z 2 ) 3/2 

where z = z (t) is the distance from the infinitesimal mass point to the plane of 
the primaries and r (t) is half the distance of the primaries, 

. 1 — ecosi? (f) 
r(*) = 2 • 

where the eccentric anomaly E (t) is the solution of the Kepler equation 

E-t + esinE, (21) 

and e is the eccentricity of the ellipses described by the primaries. We will 
assume < e < 1 all though Subsections 13.21 and 13.31 The Hamiltonian linked 
to the system is 

H = f_ 1 (22) 

2 (z 2 {t)+r*(t)) 1/2 ' 

provided v stands for z in the corresponding equations. 

Since (|22[) cannot be solved in explicit form, attempts at a Hamiltonian 
formulation of (|20[) . whether exact or approximate, require one of at least two 
options: looking for an exact Hamiltonian formulation by means of a change of 
variables, which we will do in the next paragraph, and searching an approximate 
Hamiltonian formulation, which will be done in Subsection 13.31 

Let us now find a Hamiltonian linked to [501 We may express r (t) as 

r(t) = (Ro<p)(t):= 



2(1 + ecos ¥?(£))' 
where <p, the true anomaly, is a solution of 

dip (1 + ecos f) 2 y/l — e 2 
~dt = (l-e 2 f /2 = 4i? 2 {if) ' 

and we may follow the procedure introduced in (see also [13]) by taking <f 
as the new independent variable and the new dependent variable. 

Writing t once again to denote ip, we have the following differential equation: 



ecost + (I + x 2 ) 3/2 
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clearly amenable to the hypotheses in Theorem 12.11 and in the first paragraph 
of Section [2] Defining q\ = x, q 2 = t and p\ = x, the autonomous Hamiltonian 
system corresponding to (|23[) is given by 

p? eg?coso 9 - 4 (1 + 4o?) / s 

H e = H e +p 2 :=?± + ^ ^— i ^ +p 2 , 24 

2 2(l + ecos<72j 

always assuming e £ [0,1]. 

The circular Sitnikov problem i/o is meromorphically integrable in the sense 
of Liouville- Arnold and can be solved using elliptic integrals. The non-integrability 
for e = 1 was first studied by means of straight Morales-Ramis theory in [3H 
§5] (see also [20j §5.3]); we will now extend the proof of meromorphic non- 
integrability therein to one for every < e < 1 by using Theorem 12.11 

The NVEp derived from the Hamiltonian system given by (|2"4l is 

- / e(4^ + l) 5/2 cos g2 -64^ + 8 \ 
V (eco S <z 2 + l)(4g 2 + l) 5/2 ) 

Taking q± = we have a solution T = |z (t) = ^0, t, 0, 1+e 2 cos t J | of along 
which NVEr is given by 

\ e cos t + 1 ) 

which is algebrizable, through the change of variable r = cost, into 

tL - (^-\ d ± ( eT + 8 \ £ = (25) 
dT 2 \l- T 2jdT \(er + - t 2 ) J ^ ■ V ' 

This equation can be transformed into the differential equation 

d 2 ( ( 5er 3 + 33r 2 - 2er - 30 



dr 2 V (er + 1)4(1 -t 2 ) 5 



C, (26) 



by means of £ = ^= 2 . Equations (f2"5")) and (f2"6")l are integrable in terms of 
Liouvillian solutions only if e = 0, whereas for e ^ their solutions are given 
in terms of non-integrable Heun functions if e ^ 1 and non-integrable Hyper- 
geometric functions if e = 1, i.e. their Galois groups are virtually non-solvable, 
hence virtually non-abelian; furthermore, they have a regular singularity at in- 
finity, which by Theorem 11.31 implies the non-integrability of for e S (0, 1] 
by means of meromorphic first integrals. In particular, the Galois group of 
equation (|2T))) is exactly SL2 (C). 



3.2.1 Numerical results for the Sitnikov Problem 

The author is indebted to Sergi Simon in what concerns the following subsection, 
including the figures shown at the end of the paper. Acknowledgments are also 
due to Carles Simo for further specific suggestions. 
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Let S = {sin (92) = 0}. The six figures at the end of this paper show 
Poincare sections of the flow with respect to £, projected on the (qi,pi) plane, 
for the Hamiltonian system Xh b obtained from (f!M]) . As may be easily deduced 
from said Hamiltonian, all sections are symmetrical with respect to the qi and 
Pi axes. Different amounts of initial conditions are used for the sake of clarity. 

Figure [1] corresponds to e = 0. In keeping with what was said after (f2"4")) . 
the whole subset of £ transversal to the flow sheds concentric tori (ostensibly, 
the intersections of the invariant Liouville- Arnold tori with £), a typical sign of 
integrability; the tori shown are only a selection of those therein, as the actual 
area foliated by them is larger. 

A number of these invariant tori break down upon the slightest increase in e, 
and in the ensuing figures the two most interesting features are those invariant 
sets (usually called KAM tori) whose intersection with £ prevails in the form of 
Jordan curves, and the zones of chaotic behavior between them. Sparse zones 
of the section will account for chaotic zones as well, for e > 0. Figures [2] and 
[3] show two different close-up views for the Poincare section corresponding to 
e = 0.01. The latter figure is actually a detail of the "island" of tori appearing 
at the right of the general section. For e = 0.1, Figures [4] and [5] are, respectively, 
a general view of the section and a close-up of one of the islands appearing at 
each side of the central area. As for e = 0.4, Figure [6] is an enlarged view of one 
of the two islands appearing at each side of a central area. 



3.3 The approximate Sitnikov problem 

As said in Subsection 13. 2\ we now consider an approximation of the Sitnikov 
problem; see |13] and [15j for more details. As opposed to meromorphic non- 
integrability, we will prove non- integrability by means of rational first integrals. 

The fact that <p (t) = t + O (e) yields r (t) = 1 ~ e 2 cosi + O (e 2 ) , and thus the 
Hamiltonian in (l22l) becomes 



1 cos t , 2 

H=^-^= -e ^ + 0(e 2 



whenever e ~ 0. In particular, the first-order approximation of this asymptotic 
expansion in e yields the Hamiltonian 

H^-^—^^^TH (27) 

Considering q± = x, q% = t and p\ = v; the autonomous Hamiltonian system 
corresponding to this equation is given by 



corresponding to the Hamiltonian in Theorem 12.11 

., . 8x 24a: cost 

j(x, t) = =- — e =-. 

(4x 2 + l)i (4x 2 + l)t 
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The NVE r for the Hamiltonian |2SJ is given by 

■■ = ( 24(16 g 2 -l)cosg 2 8(8g?-l) \ ^ 
V (4g? + l)i (4<7? + l)V 

Its general solution may be expressed as: 

£ (t) = K X C ^32, 48e, + ^32, 48e, | 

where the Mathieu even (resp. odd) function C(a,q,t) (resp. S(a,q,t)) is 
defined as the even (resp. odd) solution to y + (a — 2qcos (2t)) y = ([21 Ch. 
20]). 

Taking gi(t) = 0we can see that z(t) = (0, t, 0, 2e cost + 2); hence, defining 
z(t) = (0, t, 0, 2e cost + 2) and T — {z(t)}, the operator linked to NVEr is given 

by 

|= (-24ecosi-8)£, 

which is algebrizable (see Theorem ll.2j) through the change of variables r = cos £ 
into 

0-(i^)§+(^=»- 

Now, this equation can be transformed in the differential equation 

d 2 C ^96er 3 + 31r 2 - 96er - 34\ C ,, n x 



dr 2 V 4(1 — t 2 ) 2 / ^1 

Equations (f!2"9")) and (f3T)]) are integrable in terms of Liouvillian solutions only 
if e = 0, since for e ^ their solutions are given in terms of non-integrable 
Mathieu functions, hence their Galois group are virtually non-solvable and thus 
virtually non-abelian; furthermore, they have an irregular singularity at infinity, 
implying rational non-integrability for the Hamiltonian field Xg with a = 1 in 
virtue of Theorem ll.3l In particular, the Galois group of equation (f3T)]) is exactly 
SLa(C). 

Equations (j!?9"|) and (|3T)]) have been deeply analyzed in [5] using the Hamil- 
tonian change of variables t — e lt , obtaining the same result presented here. 

3.4 Painleve II equation: x = 2x 3 + tx + a 

Defining q± = x, q 2 = t, p\ = y and a £ C, the autonomous Hamiltonian system 
corresponding to this equation can given by any of the following three functions: 

H = H+p 2 , H = £-£-< b &-< m , (31) 
H = H + p 2 , J ff = ^-I( 9 2 + |) 2 - agi , (32) 
H = H+ P2 , H = &-(fi + £)p 1 -(a + ^\ 9 L 1 (33) 
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where the equations (|3"Tj) , ([3"!?)) and (|33[) correspond to the Hamiltonian of The- 
oremlO (/(x, i) = 2x 3 + tx + a), CorollaryEJ (g(x) = x 2 and a(i) = t/2) and 
Theorem 12.41 (q(x) = x 2 and a(t) = t/2) respectively. The Hamiltonian system 
for Painleve II, studied in [2U [27] , corresponds precisely to Hamiltonian ([33]) . 
The NVEr for these Hamiltonians is given by 

£ = (Qqi + ? 2 ) e 

Taking a = and gi(t) = we have particular solutions z(t) = (0, i, 0, 0), 
z(t) = (0,t,0,i 2 /8) and z(t) = (0, t, i/2, t 2 /8), respectively, for the Hamilto- 
nians ([3"Tj) . ([3^)1 and (|3"3"jl : hence NVEr is given by £ = f£, the so-called Airy 
equation ([2j §10.4.1]), which has an irregular singularity at infinity and is not 
integrable through Liouvillian solutions, i.e. its Galois group is SL2 (C), not 
virtually abelian; thus, by Theorem 1 1.3[ the Hamiltonian field Xs with a — 
is not integrable through rational first integrals. 

Now, for a = 1 and qi(t) = —l/t, the integral curve z(t) is given by 

11 1 \ /ll 1 i 2 \ , / 1 2 t 1 t 2 
respectively for the Hamiltonians (|3ip , ([32]) and ([33]) , so that NVEr is given by 

£=(& + *) t> r = {*(t)}, 



whose general solution is 



2t 3 / 2 \ (It 3 1 2 

KiI- 5/3 I — — + K 2 I 5/3 ' 



I a = 2 Q i Q ^p-pq^j + 2 2 r(2+a) ^ (^ 4 )J being, for each a, the modified Bessel 
function of the first kind, i.e. the solution to t 2 y + ty — (t 2 + a 2 ) y = ([2j 
§9.6]). 

The normal variational equation has an irregular singularity at infinity and 
is not integrable through Liouvillian functions because its solutions are given 
in term of non- integrable Bessel functions (see [23 [2?]), i-e. its Galois group 
is SL 2 (C) which is not virtually abelian; again by the Morales-Ramis Theorem 
11.31 the Hamiltonian field Xq with a = 1 is not integrable through rational first 
integrals. 
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3.5 The differential equation: x = — \ + a 

Considering q\ = x, q% = t, pi = y and a G C; the autonomous Hamiltonian 
systems corresponding to this equation are given by 

H = H+ P2 , ff=f- -f-aqi, (34) 
Z sq 1 qi 

H = H + p 2 , # = ^-1^1+2^ -aqx, (35) 

H = H + P2 , H = ^+(J-+2q 2 ^p l -(a + 2)q 1 . (36) 



(j34|) . (j35|) and ([36j) correspond to Theorem O (/(a;,t) = - + a), 

Corollary l2~3l (ofx) = -i and a(<) = -2t) and Theorem l2~4l (a (x) = -i and 
a(t) = — 2i) respectively. The NVEr for all three is given by 



3_ + 2q2 
iqf qf 



Now, for a — 1 and qi(t) — \/t, the integral curve z(t) is given by 
^=,2Vtj, (^/i,t,^=,2t 2 + 2Vi^j and (Vi,t,-2t,2i 
respectively for the Hamiltonians (|3~4"|) , ([55)1 and (|3"6")) , rendering NVEr equal to 



having a solution 



3i 3 / 2 / 7 8t 3 / 2 
?i — ^ — o-Fi 



2 7 105 v ; ' 

Fi (;a;*) = lim^oo 1F1 (q;a; |J = (^£^1 bein S tne confluent hyperge- 

ometric limit function ([H Ch. 13]), and an independent new solution £2 = 
£i/£r 2 , satisfying 

6 3^ 9 27 V ; 

As is the case for the rest of normal variational operators appearing in this paper, 
our knowledge of the exponents around of a fundamental set of solutions (in 
this case, £1 and £2), coupled with the basic result on factorization obtained 
in O Th. 8 (Ch. 5)] (see also [5J Criterion 1]) would suffice to prove non- 
intcgrability. Here, however, we will keep our restriction to Theorems 12.11 and 
and Corollarv l2.3l 
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(|37[) is algebrizable (Theorem II. 2\i , through the change of variables t = \fi 
into 

now, this equation can be transformed in the differential equation 
rf 2 C /32r 3 + 15\ A , , ,- 

Equations Ij38|) and (|39j) have an irregular singularity at t = oo and are not 
integrable through Liouvillian solutions due to the presence of Bessel functions, 
i.e. their Galois group are virtually non-solvable, therefore virtually non-abelian, 
Theorem 11.31 once again settling rational non-integrability for a = 1. In partic- 
ular, the Galois group of equation (|39j) is exactly SL 2 (C). 



4 Final Remarks: Open Questions and Future 
Work 

This paper is the starting point of a project in which the author is involved. 
The following questions arose during our work: 

• In [57] it was proven that the autonomous Hamiltonian system related 
to Painleve II is non-integrable for every a 6 Z. Is this also true for 
equation (fl3|) ? 

• Does the integrability of equation (| 13[) for arbitrary a€Z depend on the 
choice of g(x) and a(t)1 

• Assuming the above question has an affirmative answer, in what manner 
can the choice and form of g{x) and a(t) assure non-integrability for every 
q£Z? and for every a G C? 

• Is it possible to find transversal sections of the flow, and thus Poincare 
maps, for cither H or the algebraized equation, even in the absence of 
non-trivial numerical monodromies? Do Stokes multipliers contribute to 
the answer in a significant manner? 

Among our next goals, the analysis of the following items is due further 
immediate research: 

• the application of Morales-Ramis theory to higher variational equations 
of NAHS; 

• differential equations in the form x = f(x,x,t); 

• the rest of Painleve equations; 

• the theoretical aspects of NAHS such as their geometry and the feasibility 
of an analogue to Liouville- Arnold theory; 
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• the non-integrability of NAHS with two and a half degrees of freedom; 

• specific examples of NAHS related to control theory, as well as other 
related to Celestial Mechanics, such as Restricted Three- and Four-Body 
Problems and Henon-Heiles systems ([4], [11]). 

• the exact relation, perhaps causal, between separatrix splitting ([[22]) and 
non-integrability, whether rational or meromorphic. 
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Figure 2: Poincare section sin 92 = 0, e = 0.01 (Figure: Sergi Simon) 
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Figure 4: Poincare section sing 2 = 0, e = 0.1 (Figure: Sergi Simon) 
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Figure 5: Poincare section for e = 0.1 (Figure: Sergi Simon) 
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Figure 6: Poincare section, e = 0.4 (Figure: Sergi Simon) 



